Coherence versus reliability of stochastic oscillators with delayed feedback 
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For noisy self- sustained oscillators, both reliability, stability of a response to a noisy driving, 
and coherence understood in the sense of constancy of oscillation frequency belong to the main 
characteristics. Though the both characteristics and techniques for controlling them received great 
attention of researchers, owing to their importance for neurons, lasers, clocks, electric generators, 
etc., these characteristics were previously considered separately. In this paper, strong quantitative 
relation between coherence and reliability is revealed for a limit cycle oscillator subject to a weak 
noisy driving and a linear delayed feedback, a convectional control tool. Analytical findings are 
verified and enriched with a numerical simulation for the Van der Pol-Duffing oscillator. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.45.Xt 



Recently, robustness of response of a limit cycle oscil- 
lator to a noisy driving have attracted considerable at- 
tention of both experimentalists and theoreticians [l|, 0, 
0, S 0, 0, 0, !, 0, m ■ In different fields of science, related 
phenomena appear under different names. In neurophys- 
iology the reliability property of spiking neurons, which 
manifests itself as a coincidence of responses of a single 
neuron to a repeated noisy input of a prerecorded wave- 
form, attracts great attention [ij. In recent experiments 
with a noise-driven Nd:YAG (neodymium-doped yttrium 
aluminum garnet) laser a similar property has been 
referred as consistency. From the theoretical viewpoint, 
reliability and consistency are manifestations of the syn- 
chronization of uncoupled nonlinear oscillators receiving 
identical noisy driving [3, i, 0, 0, 0, 0, 0, El • 

Quantitatively, stability of response, reliability, is char- 
acterized by the largest Lyapunov exponent (LE). For 
smooth limit cycle oscillators LE is negative [1, 0, 0] 
meaning that the system is reliable. However, a large 
noise may lead to a positive LE (d, 0, 0, [n|; and even 
antireliability for neuron-like systems in a "classic" ex- 
perimental set-up has been forecasted [§]). 

However, for some oscillatory systems not only the re- 
sponse stability is important, but also the coherence, i.e., 
the constancy of the oscillation frequency, which is mea- 
sured by the diffusion constant of the oscillation phase. 
The coherence determines the precision of clocks (includ- 
ing biological ones [HI), the quality of electric generators, 
susceptibility of an oscillatory system to external driv- 
ing [13j . and predisposition to synchronization; a laser 
radiation should be coherent when one needs to focus 
the beam or redirect it without angular divergence; etc. 
In Ref. [HI (followed by methodologically closely related 
Ref. [Til) the extremely efficient technique for controlling 
the coherence by a weak delayed feedback has been pro- 
posed and theoretically analyzed (a successful experimen- 
tal implementation of this technique for a laser in chaotic 
regime has been reported in Ref. [15|). Remarkably, due 
to the time-shift symmetry a noiseless limit cycle system 
is neutrally stable and remains such in the presence of a 
delayed feedback. But in the presence of noise, the de- 
layed feedback utilized for controlling the coherence may 
considerably affect the response stability. 



Noteworthy, in the presence of both delay and noise 
(or irregularities), the process is not Markovian any- 
more; therefore, one may not apply such well-elaborated 
tools as the conventional Fokker-Planck equation, and 
ad hoc statistical methods are employed for studies 
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This paper presents both analytical 
and numerical results on the reliability of noisy-driven 
limit cycle oscillators subject to delayed feedback control, 
suggesting an effective mean for controlling the reliabil- 
ity. Analysis of these results in the context of controlling 
coherence reveals strong quantitative relations [Eq. [[15])] 
between the reliability and the coherence. The disclosed 
fact, that a high reliability occurs for a weak coherence, 
and vice versa the weaker reliability the higher coherence, 
imposes important limitations on implementation of this 
conventional control technique. Imperfect cases are also 
discussed. 

In order to demonstrate numerically the relationship 
between coherence and reliability, a simulation for noisy 
Van der Pol oscillator 

x - - x 2 )x + x = k[x(t - t) - x(t)} + e £(t) , (1) 

(€(*i)C(*a))=2<y(ti-t a ) ) (0 = o 

has been performed. Here /i describes closeness to the 
Hopf bifurcation point, k and r are the feedback strength 
and delay time, respectively, e is the noise amplitude, 
£(i) is normalized white Gaussian noise. In the presence 
of noise the oscillation phase ip = — arctan(a;/x) diffuses 
according to ((<p(t) — (</?(t))) 2 ) oc Dt. Diffusion constant 
D quantifies the coherence of oscillations. 

Fig. Q] shows the effects of a linear delayed feedback on 
the diffusion constant (DC) and the Lyapunov exponent 
(LE) measuring exponential growth rate of perturbations 
in the system ([1]). Noteworthy, not merely the LE and the 
DC are crucially magnified or suppressed simultaneously 
at t/Tq (here To is the oscillation period of the control- 
free noiseless system) being integers and half-integers, 
but even their ratio remains nearly constant as r changes 
(see Fig.Hb). 

Let us develop a phase description of the system. One 
can parameterize the states of a limit cycle system on 
the limit cycle by the oscillation phase if uniformly grow- 
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FIG. 1: (Color online) Dependencies of Lyapunov exponent 
(A) (a: upper graphs) and diffusion constant D (a: lower 
graphs) on delay time r for the Van der Pol oscillator JTJ with 
fj, = 0.7 subject to white Gaussian noise of strength e 2 = 0.01 
and the linear delayed feedback of strength k — 0.06 (squares) 
and k = —0.06 (circles). Oscillation period of the control- free 
noiseless system To ~ 27r/0.96 . The solid lines present ana- 
lytical dependencies [Eqs. Qli p .l|12 p ]. (b): The inconstancy of 
ratio ~{\)/D is not resolvable against the background of the 
calculation inaccuracy. 



ing in the course of temporal evolution. Such an oscil- 
lator subject to weak noise and feedback stays in the 
vicinity of this cycle, and its evolution may be still de- 
scribed within the framework of the conventional phase 
approximation [20]. Close to the bifurcation point, i.e., 
for /i — > 0, the Van der Pol oscillator has a nearly circular 
limit cycle: xq = 2 cos ip, Xq = —2 shop, and the phase 
equation for the system (TT]) reads (cf. 



(p = n + ag[<p(t - r), <p(t)] + ef[<p(t)] ° , (2) 

where fio = 27r/To is the inherent cyclic frequency of the 
system, a — fc/2, g — sin[</?(t — t) — <p(t)], the sign "o" 
means a Stratonovich form of the equation, 27r-periodic 
function /(</?) is the sensitivity to noise. For an additive 
noise as in Eq. {!]), f(ip) = (2fi ) — 1 cos<£> (c f. I13ll ). but 
we keep / for generality. Noteworthy, in Ref. [18[ , Eq. @ 
has been used to describe the evolution of the phase of 
an optical field in a laser with a weak optical feedback. 
For a small perturbation a, one finds 

a = acos[<p(t — t) — ip(t)](a(t — r) — a(tj) 

+ ef'[<p{t)]a{t)offl 

(the prime stands for derivative with respect to the ar- 
gument). Therefore, the instant exponential growth rate 
X(t) obeys 

A(t) = a cos[ip(t - t) - tp(t)] (e - A(tl)dtl - l) 

+ e/V(t)W(t), (3) 
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FIG. 2: (Color online) Same dependencies as in Fig.fj] for 
the same parameter values but for red Gaussian noise £(t) = 
T~ J\_ T £(ii) dti with T = 1.2 . For notation see caption to 
Fig-IB 



here we have made use of a(t) oc exp(J* A(ti)dti). Note, 
that the LE is the mean value (A). 

For further analysis, it is more convenient to consider 
equations in Ito form. Eqs. (HI),© read 

= Cl + aam[<p(t-T)-<p(t)]+ e 2 f f+ ef (4) 

A(i) = a cos[ip(t - t) - <p(t)] (e~ K-* A(tl) dtl - l) 

+ e 2 /"/ +£/'[¥>(*)]£(*)• (5) 

The terms ahead of the noisy ones describe the 
Stratonovich drift. Recall, in Ito form (with the 
Stratonovich drift included explicitly) the instant value 
ip(t) is independent of the instant value taken at the 
same time moment t. 

Let us explicitly introduce the mean frequency O and 
the instant frequency deviation v; ip = fit + ip, if) = v, 
(v) = 0. For a weak noise and a small feedback strength 
(e <C 1, \a\ <C 1) the instant frequency fluctuations are 
small (v <C 1), and Eqs. (HI), (HI yield, up to the main 
order of accuracy, 

fl = Oo — a sin fir , (6) 
v = -a cos nr [^(t- T )-^(t)} + e 2 f f + ef £(i), (7) 

A = a(cosi7r + sin fir [ip(t — r) — ip(t)f\ 

x (e" A-r dt i - l) + £ 2 f"f + ef £(t) . (8) 

Let us now find the LE. Assuming v and the fluctuating 
part A of A, obeying 

X(t) w -acosOr f A(*i) dt x + ef'(ttt) ^(t) , (9) 

Jt-T 

to be Gaussian, one can employ Furutsu-Novikov for- 
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FIG. 3: (Color online) Dependencies (A)(r) (a: upper graph) 
and D(t) (a: lower graph) for the Van der Pol-Duffing os- 
cillator dHJ) with fj, = 0.2 (T ft* 2tt/2.02) subject to white 
Gaussian noise, e = 0.05, and the delayed feedback, k = 0.06. 
For description and notation of (b), see Fig. [Tj caption. 



mula [HI to obtain from Eq. © |23( 

(A) = acosr!r(- (A)r+i([| t A(ti)d*i] 2 )) ~e 2 (f' 2 ) v 

(10) 

(here (•••)^ stands for an average over the phase if). 
The value I = / [//_^(ti) d*i] 2 ) can be evaluated from 
Eq. © (similarly to (u(ii)u(* 2 )) in Ref. [13]); 



/ = £ 2 <f' 2 },- 



£ 2 (/' 2 ), 




IX 



1-e" 
2r 



1 + ar cos fir 



1 + ar cos Or 



- cos Or dx 
for ar cos Or > — 1 ; 
for ar cos Or < — 1 . 



For ar > 1, Eq. © exhibits multistability of mean fre- 
quency O [lT], [l8|], which results in the violation of 
the basic assumptions of our analytical theory. Hence, 
the case ar cos Or < —1 may be ignored as meaningless, 
and after substitution of /, Eq. (|10p reads 



(A) 



e 2 (f% 



(11) 



(1 + ar cos Or) 2 ' 
whereas the DC has been already evaluated in Ref. [HI]; 



D 



2e 2 (/ 2 }, 



Therefore, 
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FIG. 4: (Color online) Same dependencies as in Fig. [3k for 
the same Van der Pol-Duffing system and feedback strength 
but a larger noise, e = 0.374 . 



which is 1/2 for f(tp) = (20o) _1 coscp as for the Van der 
Pol system in Fig.[T] Note that, due to the deformation of 
the limit cycle at /i = 0.7, relation (fl3|) is more accurate 
than Eqs. (fTTj) . (fT2|) where the term ar cos Or is specific to 
g = sin[tp(t-T)-<p(t)] (see Fig.[[]where -(X)/D w 0.55). 

While deriving Eqs. (fTTj) , (fl~2"j) we nowhere utilized that 
the noise is (^-correlated. Remarkably, the results remain 
valid for colored noise, e.g., a red one (see Fig. [J), re- 
sults for which coincide with the one for white noise 
almost up to the numerical calculation inaccuracy and 
-(A red )/£> re d « 0.53. 

For a strong noise the phase description always leading 
to a negative LE is not applicable, and even positive LEs 
have been reported [E B H B • This case can be treated 
only numerically. For this sake, a simulation for the Van 
der Pol-Duffing oscillator 

x- n(l-x 2 )x + x + x 3 = k[x(t-T)-x(t)]+e((t) (14) 

exhibiting positive LEs for a moderate noise Q has been 
performed. Let us note that for a non-large noise (Fig©, 
ratio — (A)/£> is changed not greater than by 20%, while 
the DC and the LE are changed by factor w 20, in a 
broad range of r (the only exception is interval [To/2, To] 
near the domain where the LE is positive). 

At e — 0.374 (Fig©, where the control- free Van der 
Pol-Duffing oscillator just becomes unstable (unreliable), 




FIG. 5: The snapshots of the ensemble of 20 000 identical 
Van der Pol-Duffing oscillators (/i = 0.2) driven by common 
white Gaussian noise (e = 0.374) have a fine structure for r 
indicated in the plots (cf. to the LE in Fig. [3}. 
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(A) \ T= o « 0, the linear delayed feedback leads to maximal 
positive LEs for integer t/Tq and minimal (not always 
negative) LEs for half-integer t/Tq. Concerning the in- 
terpretation of dependence (A) (r) , note the following. In 
the absence of the feedback control, an intermittency of 
epochs of positive and negative local LEs ( "local" means 
evaluated over a finite time interval) takes place and the 
transition to positive LE is related to a plain quantita- 
tive prevalence of the former over the latter (cf. [9| ) . The 
feedback affects (magnifies or suppresses) the local LEs 
over these epochs non-uniformly, thus shifting balance 
between these epoches and bringing about a domination 
of positive local LEs for integer t/Tq and negative ones 
for half-integer t/Tq. For positive "global" LEs, phase 
diffusion owes mainly not to stochasticity but to chaos 
(samples of snapshot chaotic attractors, see Ref. (2^, are 
presented in Fig. [5]). As a result, here the DC is dimin- 
ished where the LE is minimal. 

Summarizing, for a weak white or colored Gaussian 
noise [2J], highly stable response (reliability) to a noisy 
driving is observed when phase diffusion is strong (i.e., 
the coherence is weak). Vice versa, for small diffusion 
(i.e., highly coherent oscillations) response is weakly sta- 
ble [Figs. 1,2, Eq. (fl3|) ]. In particular, this imposes strong 
limitations on the implementation of the technique of co- 
herence improvement by virtue of a linear delayed feed- 



back. For instance, in an ensemble of uncoupled iden- 
tical self-sustained oscillators synchronized by a com- 
mon external noisy driving, small intrinsic noise is al- 
ways present and leads to spreading of oscillator phases: 
Aip oc ejn/ yj — (A) (e; n is the amplitude of intrinsic noise, 
cf. [6[). In such an ensemble the delayed feedback im- 
provement of the coherence results in a mutual spreading 
of oscillator phases which may be sometimes undesirable. 

For a strong noise being capable to create a positive 
Lyapunov exponent, i.e., antireliability, chaotic contri- 
bution to phase diffusion may prevail over the stochastic 
one, and then an enhanced coherence occurs for the max- 
imal reliability (Fig.|4]). 

Detailed calculation of the Lyapunov exponent in the 
above text serves the purpose to disclose the essentially 
different nature of various contributions to the Lyapunov 
exponent and the diffusion constant. However, the final 
quantitative effect of delayed feedback on these dissimilar 
properties of oscillatory systems somewhat surprisingly 
turns out to be identical. The reported phenomenon be- 
ing valid for a general class of limit cycle oscillators is, 
thus, neither intuitively expected nor trivial. 
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